Opinion dynamics on directed small-world networks 
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In this paper, we investigate the self-aSirmation effect on formation of public opinion in a directed 
small- world social network. The system presents a non-equilibrium phase transition from a consensus 
state to a disordered state with coexistence of opinions. The dynamical behaviors are very sensitive 
to the density of long-range interactions and the strength of self-affirmation. When the long-range 
interactions are sparse and individual generally does not insist on his/her opinion, the system will 
display a continuous phase transition, in the opposite case with high self-affirmation strength and 
dense long-range interactions, the system does not display a phase transition. Between those two 
extreme cases, the system undergoes a discontinuous phase transition. 
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I. INTRODUCTION 



II. MODEL 
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Recently, much effort has been devoted to the studies 
of opinion dynamics [il,J3, 0. The Sznajd model [ii], 
Galam's majority rule [g, |7| , and the Axelrod multicul- 
tural model [s] describe opinion dynamics as individu- 
als follow their nearest neighbors' opinion. However, the 
real-life system often seems a black box to us: the out- 
come can be observed, but the hidden mechanism may 
be not visible. If we see many individuals hold the same 
opinion, we say The Spiral of Silence phenomenon oc- 
curs. It is common in real world that people adhere to 
their own opinion even opposite to most of their friends 
[Tol . [Til . [T3 |. which we call self-affirmaUon of individu- 
als. Traditional opinion models fail to account for the 
self-affirmation effect of individuals. In our early work, 
the influence of inflexible units has been investigated in 
a simple social model 13]. It is found that this kind of 
effect can lead to a nontrivial phase diagram. 



In the real world, interactions between individuals 
are not only short ranged, but also long ranged [l^, 
[Tsj . The interaction usually takes on directed feature, 
which means an individual who receives influence from 
a provider may not effect the provider. We use directed 
links to represent this kind of relation between individ- 
uals in the directed small-world networks proposed by 
Sanchez et al pj6j . In this paper, we present an opinion 
dynamics model including individual self-affirmation psy- 
chological feature and directed long-range correlation be- 
tween individuals. The parameter space can be roughly 
divided into three regions, in which, respectively, we ob- 
served continuous phase transition, discontinuous phase 
transition and no phase transition. 
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In this section, we introduce a directed small-world 
network model and an opinion dynamics model. We start 
with a two-dimensional regular lattice, in which every 
node is connected with adjacent four nodes inwardly and 
outwardly respectively, then, with probability p, rewire 
each link connected outwardly with a neighbor node to 
a randomly chosen nonadjacent node. In this way, as 
shown in Fig. 1, a directed network with a density p of 
long-range links is obtained. In the network, nodes rep- 
resent individuals in the social system and the outward 
links represent the influence from others. Each node con- 
nects with four nodes outwardly which are called mates 
of the node. Long-range links are used to describe the 
long-range correlations between individuals. 

In the network, state of each node represents the view- 
point of the corresponding individual, which evolves ac- 
cording to social process, determined not only by other 
correlative surrounding effects but also by its own char- 
acter. To illustrate these cases simply, it is supposed 
that there are two kinds of possible opinions in the 
system, just as the agreement and disagreement in the 
election, and each individual takes only one of them. 
Therefore, the state of a node i can be described as ct^, 
CTi G {+1, —1}- We describe the difference of (Ji between 
its mates' state by the function W{ai) ^ 2ai J2'j=i '^j' 
which (Tj (j = 1, 2, 3, 4) are states of mates of node i. In 
addition, q (0 < g < 1) is used to describe the probabil- 
ity, with which individuals follow their mates' dominant 
opinion. Meanwhile I — q represents the self-affirmation 
probability of individuals, with which an individual in- 
sists on his/her own opinion though it is opposite to most 
of his/her mates. 

According to the illumination above, we introduce the 
dynamic rule as follows: W{cri) > indicates that dj is 
the same as the majority of aj {j = 1, 2, 3, 4) and ai over- 
turns with probability exp[—W{ai)/T] which depends on 
a temperature-like parameter T. W{ai) < represents 
that ai is opposite to the majority of aj (j = 1,2,3,4), 



FIG. 1: Illustration of the structure of a directed small- world network for p = 0.1 [l6l ). 



and ai overturns with probability q. When W{ai) = 0, 
we consider that the state of node i overturns with prob- 
ability q also. So that the overturning probability P{<Ji) 
of Gi is given by 

r exp[-M^(a,)/r] ,forW^(a,)>0 
P^^^)-]^ q , for < 0. ^ 

From the dynamical rule (1), we can see that, when q = 1, 
the current model restores to the network-based Ising 
model However, our model is non-equilibrium be- 

cause the overturning probability of a state does not sat- 
isfy the detailed equilibrium condition. 



III. SIMULATIONS 

In order to describe the evolution process of the model, 
we employ a magnetization like order parameter 



™=l e {+1,-1}. (2) 

1=1 

The network size is L x L and m represents the absolute 
average value of the states of all nodes. An extensive 
Monte Carlo numerical simulation has been performed 
on our model with a random initial configuration and a 
periodic boundary. Results are calculated after the sys- 
tem reaches a non-equilibrium stationary state. In order 
to reduce the occasional errors, for network size L = 16, 
32, 64, and 100, we have averaged the result over 40000, 
10000, 2000, and 1000 runs, respectively, with different 
network structures under different random initial config- 
uration. Obviously, when (m) tends to 1, the system en- 
ters into an ordered state, i.e. individuals in the system 
reach the opinion consensus. Meanwhile, if the system 
stays in a disordered state, the order parameter scales as 
(to) ~ As shown in Fig. 2, the system reaches an 
ordered state when T is less than a critical temperature 
Tc- The system displays a continuous phase transition 
for p — 0.1 and q = 0.9, while a discontinuous phase 
transition for p = 0.9 and q — 0.9. From the probability 
density functions (PDFs) of the order parameter near the 



phase transition point of the phase diagram (p, T) , one 
can distinguish between the continuous phase transition 
and discontinuous phase transition clearly. According to 
PDFs inserted in the upper and lower of Fig. 2, It is 
found that the most probable values of to, which cor- 
respond to the highest peaks of PDFs, jump little from 
nonzero to zero in continuous phase transitions from Fig. 
2a to Fig. 2b, while largely in discontinuous ones from 
Fig. 2c to Fig. 2d. It seems that the long-range correla- 
tions can change the nature of phase transition. 

Evidently, the system varies from the continuous phase 
transition to discontinuous phase transition when the 
density of directed long-range connections is high enough 
for q = 0.9. It is natural to ask how these topology 
structures influence the opinion dynamics. To solve this 
problem, we define the domain size s as the number of 
neighborhood nodes in the same state. As showed in 
Fig. 3a, it is found that the domain size s distributes in 
a power law, g{s)'^s~'^ for s <IC at the critical point, 
where g{s) is the probability function. One can find that 
there is a local maximum probability of large domain for 
p = 0.1 compared to the size distribution for p = 0.9 in 
Fig. 3a. Smaller p indicates more localized interactions 
between individuals, and a large domain emerges more 
easily. Besides, we calculate the number of time steps, i, 
during which an individual holds the same opinion. As 
shown in Fig. 3b, one can find that individuals change 
their own opinion for p = 0.9 more frequently than for 
p ~ 0.1 at T = 0.1, and the probability of t obeys a 
power-law distribution f{t)^t~'^ {t < to) for p — 0.1. 
Clearly, p plays the key role in determining the commu- 
nication strength between different opinion domains, and 
individuals change their own opinions more frequently 
because of the effect of long-range connections between 
different opinion domains. 

Fig. 4a, 4b, and 4c show the phase diagram of opinion 
dynamics determined by the network structure parame- 
ter p as well as the individual self-affirmation psychology 
characteristic parameter 1 — q. In Fig. 4a for q=0.9, the 
system displays continuous phase transition for p < pc, 
while discontinuous phase transition for p > Pc- The sys- 
tem displays discontinuous phase transition for g = 0.5 in 
Fig. 4b. The system displays discontinuous phase tran- 
sition for p < po and q = 0.3 in Fig. 4c, while the system 
does not have a phase transition for p > po and q = 0.3 in 
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FIG. 2; (m) varies with T for different system sizes. The upper and lower plots are for p = 0.1 and p — 0.9, with q = 0.9 fixed. 
Insets are PDFs nearby the phase transition point; (a) T T' , (b) T ~* T+ , (c) T ~* T;r , (d) T ^ r+. 




FIG. 3: (Color online.) The distribution of domain size g{a) (a) and opinion holding time f{t) (b) in different networks with 
q = 0.9 fixed, (a) for T = Tc and (b) for T = 0.1. The data points are obtained from 10"" samples with fixed network size, 
L = 64. 
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Fig. 4c. As shown in Fig. 4d, the continuous phase tran- 
sition takes place in the area /, the discontinuous phase 
transition appears in the area // and the system stays in 
disordered state without phase transition in the area ///. 
When both the parameters p and 1 ^ q arc large enough, 
indicating weak interaction between individuals in both 
local and global levels, the system keeps in disordered 
completely at any temperature, i.e. the phase transition 
can not take place in the system, as in the area /// of 
Fig. 4d. 

A finite-size scaling analysis is employed to study the 
critical behavior of continuous phase transition for p — 
0.1 and q = 0.9. In the neighborhood of the critical 
point Tc, (m) oc (Tc ~ Tf , (T < Tc), where /? is the 
order parameter exponent. Besides, when T is near to 
critical point Tc of the second order phase transition, a 
character length scale ^ denotes the correlation length in 
space. ^ oc (Tc — T)~'^, (T < Tc), where is a correlation 
length exponent in the space direction. At critical point, 
various ensemble-averaged quantities depend on the ratio 
of system size and the correlation length L/^. Therefore, 
the order parameter (m) satisfies the scaling law in the 
neighborhood of the critical point: (m) oc L~^^'^ f[{Tc — 
T)L^/'']. At Tc, (m) oc L-'^/", and we obtain P/u = 
0.530(5) for p = 0.1 and q = 0.9 in Fig. 5(a). Fig. 
5(b) reports (m)i^/'' versus (1 - T/Tjii/"^ on a double- 
logarithmic plot for g = 0.1 and q=0.9. It is shown that 
with the choices 13/v = 0.530(5) and v = 0.92(1) the 
data for different network sizes are well collapsed on a 
single master curve [l3|- The slope of the line is /3 = 
0.488 ± 0.005, which gives the asymptotic behavior for 
(m as L — > oo. So that, we have (3 = 0.488(5), 

V = 0.92(1). 

IV. CONCLUSION 

In conclusion, the effect of directed long-range links 
between individuals on the opinion formation is system- 



atically explored. The results show that the system takes 
on a non-equilibrium phase transition from a consensus 
state to a state of coexistence of different opinions. With 
increasing density of long-range links, a continuous phase 
transition changes into a discontinuous one. The reason 
why the phase transition behavior varies is that the long- 
range links make individuals change their own opinion 
more frequently. It is worth mentioning that the sys- 
tem keeps in a disordered state when there are sufficient 
long-range links. Those long-range interactions break the 
possibly local order, thus hinder the global consensus. 

In opinion dynamics, the self-afhrmation psychology 
character sometime may lead to polarized decision |18l . 
[T9| . Moreover, interaction between individuals in so- 
cial srotem depends on the topology of social networks 
[20I . [2H . In macroscopic level, the opinion dynamics is 
highly affected by social structure, while in the micro- 
scopic, it is sensitive to the dynamical mechanism of in- 
dividual. Our work shows a systematic picture of opinion 
dynamics, and provides a deep insight into effects of these 
two factors. 
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FIG. 5: Finite size scaling of continuous phase transition for p = 0.1 and q — 0.9. (a) A log-log plot of the order parameter 
(m) against L. (b) Double logarithmic plot of (m) L''^" versus (1 - T/T^)L'^''' for L = 16, 32, 64, and 100. 
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